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METHOD AND SYSTEM FOR AUTOMATIC WATER 
DISTRIBUTION MODEL CALIBRATION 

CROSS REFERENCE TO RELATED APPLICATION 

This application claims the benefit of the filing date of the following application 
of which it is a continuation: United States Patent Application Serial No. 09/990,81 8, 
filed November 14, 2001, which is incorporated herein by reference. 

BACKGROUND OF THE INVENTION 

FIELD OF THE INVENTION 
This invention relates to calibration of a water distribution model. 

BACKGROUND INFORMATION 

In many civil engineering applications, and ui particular, in water distribution 
networks for cities, townships and municipalities, it is of the utmost importance to have 
predictive simulations for the water pipe flow and hydraulic pressure conditions for the 
water network in that area. This is to ensure the availability of the potable water resource 
to the community, as well as to be able to perform vulnerability studies to produce risk 
assessments for risks that may be presented, which could affect the water supply. 

A water distribution model is employed for these purposes in which data describ- 
ing field conditions are assembled in a systematic manner to describe pipe flow and junc- 
tion hydraulic grade lines (HGL), which are pressures, within the water network. The 
model is desirably capable of simulating conditions encountered at the site. 
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Engineers and modelers often calibrate the models they construct. A calibration 
involves the process of adjusting model characteristics and parameters so that the model 
predicted flows and pressures match actual observed field data to some desirable or ac- 
ceptable level. Model calibration would not be necessary if the mathematical model ex- 
actly represented the actual physical processes and there was a perfect knowledge of all 
required parameters. Unfortunately, these criteria are not generally met; thus calibration 
is considered the most desirable way to achieve accurate model adjustment to most 
closely represent the water network being simulated. 

Traditionally, methods for calibrating water distribution network models rely 
upon field measurement of network pressures, pipe flows, and water levels in storage fa- 
cilities. The model is constructed, and then field observations are made by an engineer 
who visits various locations on the site and takes field observed measurements of pipe 
flows, water levels and pressures. The model is then adjusted on a trial and error basis so 
that the model simulation result more closely represents the observed data. 

For example, a model representation of a water network may be developed which 
may include information about 12-inch mains, major 8-inch mains and loops and pipes 
that connect to sampling sites. A roughness coefficient is assumed for all pipes. Another 
aspect of the network is that of "demand." "Demand " as used herein, relates to the con- 
sumer demand for water at a given point in time. Demand patterns can be estimated 
based on the number of structures of diflferent types in conjunction with an average water 
use by structures. Using this information, a working model is built to indicate how the 
network will behave in the real world to determine, for example, how much water is used 
at certain points in the network. Information is inserted into the working model, such as 
pipe roughness, and a basic working model is produced. 

At this point, the model is then calibrated. As noted, prior techniques involved a 
trial and error process by which an engineer or modeler monitors various values such as 
pressure and flow to obtain a predicted model to compare to the observed data. If the 
predicted model does not compare closely with the observed data, the engineer retums to 
the working model, makes some adjustments, and runs it again to produce a new set of 
simulation results. This may have to be repeated many times to make sure that the 
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working model produces a close enough prediction of water network behavior in the real 
world. 

There are several disadvantages to the traditional calibration methods. One such 
disadvantage is, in a steady-state simulation, it is desired to provide no changes during 
field observation in the relevant aspects of the network. And yet, the observation itself 
coxxld incorporate data from different network states. More specifically, an engineer or 
perhaps several engineers, take measurements in the field sequentially. During the time 
elapsed between taking the various measurements, the state or condition of certain as- 
pects of the network can change. A simple, but illustrative, example is that of an engi- 
neer measuring pipe flow at location A, at which time a network pump may be in an 
"ON" position, thus the pump is operating. Later, when the engineer takes a field ob- 
served measurement at location B, the pump may now be in an "OFF" state, which would 
change pipe flow (and pressure) readings within the network. The network has changed 
during the observance of the data in the field, thus affecting the accuracy of the results. 

A further disadvantage of traditional modeling techniques is that they are, among 
other things, quite time consuming. A typical network representation of a water network 
may include hundreds or thousands of links and nodes. Ideally, during a water distribu- 
tion model calibration process, the roughness coefficient and pipe diameter is adjusted for 
each link, and demand adjusted for each node. Typically, however, only a percentage of 
representative sample measurements are used in a model, due in part to the time and labor 
requirements associated with gathering the evidence. 

In addition, the model calibration process conventionally used does not take into 
account user weighted observation data such that the user can adjust hydraulic grade line 
(HGL) and/or pipe flow at data points of particular interest or importance to the user. 
Fxirthermore, known model calibration techniques use only one input parameter, pipe 
roughness, and this parameter is typically not weighted for the particular network in- 
volved. 

In addition, when model calibration software is employed the software is run and 
a set of results is produced. However, the user cannot terminate or pause the application 
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during run time to observe data at particular points, but instead, must wait until the appli- 
cation runs in full to then observe a single calibration solution so produced. 

There remains a need, therefor, for a calibration process that results in a highly 
accurate model of a water distribution network. There remains a further need for such a 
process that does not involve undue amounts of trial and error in which multiple moni- 
toring and measuring visits to the site must occur. There remains a further need for a 
method which produces a more reliable model, and allows the user to employ a number 
of weighted parameters which more accurately reflect the particular network being mod- 
eled in order to customize that model so that it more closely represents the actual behav- 
ior of the network. 

There remains a further need for a modeling system in which the calibration can 
be performed automatically and which calibration process can be refined and manually 
adjusted during the calibration run time. 

It is an object of the present invention to provide a calibration system that 
achieves these goals and that includes automatic calibration that takes into account a 
number of parameters and boundary conditions. 

SUMMARY OF THE INVENTION 

The disadvantages of prior techniques have been overcome by the present inven- 
tion, which provides an automatic water distribution model calibration process that has 
improved accuracy and efficiency. Multiple parameters and corresponding boundary 
conditions are taken into account to provide an accurate representation of the network at 
an uistant in time. The inventive system includes a software program that contains three 
integral parts: a genetic algorithm module, a hydraulic simulation module and a calibra- 
tion module. These modules interact to provide an optimized calibration solution. 

More specifically, the invention includes a method of automatically calibrating a 
water distribution model that involves a user selecting multiple calibration parameters. 
For example, parameters may be chosen including the pipe roughness coefficient, junc- 
tion demand, and pipe and valve operational status, or any combination of the parameters. 
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Next, the user enters field observed data, namely amounts for pressure and pipe flow. 
The observed data can be weighted, in accordance with the invention, with a user- 
selected weighting function to focus the calibration on certain data points. The demand 
loading information is required for different times of day corresponding to the time when 
the observed data is collected. For example, there will be a higher figure for the demand 
loading for high demand times, such as in the moming, and lower amounts for lower de- 
mand times such as in the overnight hours. Boundary conditions such as: storage tank 
levels, pressure control valve settings and pump operation speeds are also entered. In this 
manner, accuracy is improved by providing a realistic snapshot of the network actually 
operating at each instant in time. 

A genetic algorithm program then generates a population of trial solutions of the 
model calibration. A hydraulic network solver program then simulates each trial solu- 
tion. More specifically, the resulting hydraulic simulation predicts the HGL (junction 
pressures) and pipe flows at a predetermined number of nodes (or data points) in the net- 
work. This information is then passed back to the associated calibration module. 

The calibration module evaluates how closely the model simulation is to the ob- 
served data. In doing so, the caHbration evaluation computes a "goodness-of-fit" value, 
which is the discrepmicy between the observed data and the model predicted pipe flows 
and junction pressures or HGL, for each solution. This goodness-of-fit value is then as- 
signed as the "fitness" for that solution in the genetic algorithm program. One generation 
produced by the genetic algorithm is then complete. The fitness measure is taken into 
account when performing the next generation of the genetic algorithm operations. To 
find the improved "survivors," in this case, the optimal calibration solutions, fitter solu- 
tions will be selected by municking Darwin's natural selection principal of "survival of 
the fittest". The selected solutions are used to reproduce a next generation of calibration 
solutions by performing genetic operations. Over many generations, the solutions evolve, 
and the optimal or near optimal solutions ultimately emerge. 

Multiple near optimal solutions can be made available at the end of the genetic al- 
gorithm run. The best solution can be kept, and a number of top solutions produced are 
also kept. During run time of the optimization, the user can pause the program and ob- 
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serve values at that point, as desired. After the optimization step, a number of model pa- 
rameter solutions are produced^ The user can apply engineering judgment to choose one 
of these solutions to build a calibrated model, and can also then make manual adjust- 
ments to it to refine the results. 



BRIEF DESCMPTION OF THE DRAWINGS 

The invention description below refers to the accompanying drav^ngs, in which: 

Fig. 1 is a screen shot of a water distribution network to be modeled and cali- 
brated using the method and system of the present invention; 

Fig. 2 is a pipe section that may be employed in the network of Fig. 1; 

Fig. 3 is a schematic illustration of the data flow in accordance with the present 
invention; and 

Figs. 4A and 4B depict a flowchart of a procedure followed in performing the 
method of the present invention. 

DETAILED DESCRIPTION OF AN ILLUSTRATIVE 

EMBODIMENT 

Fig. 1 is an illustration of a screen shot of a user interface employed as part of the 
present invention in which a water distribution network 100 is depicted schematically. 
The water distribution network has individual main lines 104, 106, which may include 
ten-inch diameter pipes, eight-inch diameter pipes, and other hydraulic components. 
Within the network, a number of pump stations, and valves (not shown) will also be in 
operation to move and control the flow of water in the network. The pipes, valves and 
pumps are sometimes collectively referred to herein as "links." 

A pipe portion is illustrated in Fig. 2. In modeling and calibrating the model of a 
water distribution network, in accordance with the present invention, a number of pa- 
rameters are employed to develop predictions of how the network will behave in opera- 
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tion in the real world. Pipe roughness is one such parameter, such as the roughness 
schematically illustrated in pipe portion 200 as area 204. The roughness coefficient, and 
pipe diameter 206 are mathematically related to the velocity of flow 208 through the pipe 
200 using the following equation: 



where F represents the velocity of pipe flow, D is the pipe diameter, A/ is the hy- 
draulic head loss, /is the pipe friction coefficient and L designates the pipe length. The 



rameter in the modeling process. 

Other parameters include junction demand, which is the demand for water at an 
instant in time at a given junction 108 on the network (Fig. 1). As noted, other parame- 
ters include junction pressure, as well as pump and valve operational status. 

Fig. 3 is a schematic illustration of the data flow through software modules using 
the application and user interface embodying the present invention. A user interface 302 
on a personal computer or other workstation lends the user the ability to enter the field 
observed data 304, the demand alternatives 306 and other information 308 for the net- 
work. Next, the user chooses the calibration parameters such as pipe roughness, junction 
demand and link status, which data is schematically illustrated in block 310. 

The field observed data is next transferred to a calibration software module 308. 
An initial calibration model is then established. It can be either presented to a genetic 
algorithm module 3 1 1 for automatic calibration and/or to a hydraulic simulation model 
3 12 for a manual calibration run. Boxmdary condition information corresponding to the 
state of the system, at each point m time when the field observed data was collected is 
also introduced and entered into the calculation. This produces results that more accu- 
rately reflect the system being modeled. 

As noted herein, a genetic algorithm is a known model of machine learning, 
which is derived from a metaphor of the processes of evolution in nature. This is done by 
the creation (vrithin a machine) of a population of individuals represented by chromo- 




pipe friction coefficient is related to the pipe roughness, hence pipe roughness is a pa- 



7 

\\CHEETAH\VOLl\CLIENTS\107\051\0001Cl\PROSECUT\latest patapp.doc 01/17/02 9:18 AM 



PATENT 
107051-OOOlCl 

somes. In essence, it includes a set of character strings that is analogous to the base 4 
chromosomes seen in DNA. The individuals then go through a process of evolution. 
Genetic algorithms can be used for an application to obtain a multi-dimensional optimi- 
zation of a problem m which a character string (the chromosome) can be used to encode 
values for different parameters being optimized. 

In the present invention, the genetic algorithm generates a population of trial so- 
lutions of model calibration. Each trial solution is then presented to a hydraulic network 
simulation module 3 12 of the present invention. The hydraulic network simulation mod- 
ule 312 embodies a hydraulic network solver that runs a hydraulic simulation to predict 
the junction pressures and pipe flows at nodes in the network 1 00. This information, in 
the form of predictions, is passed back to the calibration module 308 of the present in- 
vention. 

The calibration module 308 then computes a "goodness-of-fit" value, which is the 
discrepancy between the model predicted and the field observed values originally input 
by the user at 302. Using the "goodness-of-fit" information, the genetic algorithm opti- 
mizer 3 1 1 reproduces better solutions, as the solutions evolve generation after generation. 
The optimal solution is found at the end of the genetic algorithm run. The fitness of all of 
the individual parameters in the population is evaluated. A new population is obtained by 
forming genetically suggested operations such as crossover, reproduction or mutation, on 
the individuals whose fitness has just been measured. The old population is discarded 
and the new one is itemted. The first generation of the process operates on the population 
of randomly generated individuals. From thereon the genetic operations in concert with 
the fitness measure operate to improve the population and select the "survivors." In the 
present application, the survivors are the calibration parameters to be used in the water 
distribution model constructed firom the information contained in database 3 14, to most 
accurately reflect the behavior of the network in the real world in its physical implemen- 
tation. 

The flowchart of Figs. 4A and 4B represent a procedure 400 representing the 
method of the present invention that begins at the Start step, designated by reference 
character 402. The user supplies the field observed data and selects parameters such as 
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pipe roughness, junction demand and link statxis as shown in step 404. As noted, the user 
may deliberately weight the observed data for focussing the calibration on critical data 
points. In order to do so, the user selects one of four weighting functions which can in- 
clude, but are not limited to, linear, square, square root and logarithm to conduct the cali- 
bration on the weighted, observed hydraulic grade line (HGL) and/or pipe flows, as 
shown in steps 406 and 408. It may be determined that other weighting functions may 
also be used while remaining v^ithin the scope of the present invention. 

The pipe roughness can be either computed as a new roughness value for a group 
of pipes within tiie network 100, or modified by multiplying the initial roughness v^th a 
correction factor. The junction demand can be selected as a calibration parameter as 
well. Adjusting the junction demand enables the user to calibrate the model to better re- 
flect the real demand condition when the observed data is collected, thus it improves the 
accuracy of the overall model calibration. Automatically identifying the link status (open 
or closed) enhances accuracy, but also assists engineers in quickly clarifying fault opera- 
tion settings for valves, pipes and pumps in order to perform diagnostics and trouble- 
shooting of the network. 

As noted, loading and boundary conditions are also supplied. Multiple loadings 
can be employed, and a number of multiple boundary conditions may also be used, as 
discussed previously herein. Now that the initial information has been entered into the 
system, the genetic algorithm module generates a population of trial solutions of model 
calibration, step 410. 

The hydraulic network simulation module then runs each trial solution, step 412. 
The hydraulic simulation predicts the HGL (junction pressures) and pipe flows at various 
nodes in the network as illustrated in step 414. Thereafter, as illustrated in step 416, the 
discrepancy between the model predicted and the field observed values of HGL and pipe 
flow are measured. The "goodness-of-fit" is computed for each solution as shown in step 
418 of Fig. 4B. The "goodness-of-fit" is assigned as the "fitness'' of that solution for 
purposes of the genetic algorithm, step 420. Using these entries as the fitness of each 
solution, the genetic algorithm operations are then again employed (in a repetitive opera- 
tion), as shown in step 422 and by the feedback loop to the step of block 410, to repro- 
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duce each next generation of calibration solutions. It is the genetic algorithm that gener- 
ates and searches for an optimal solution or several top solutions, step 426, that corre- 
sponds with the minimum discrepancy between the simulated and the observed HGL and 
flow data. 

During calibration run time, the user may terminate the program or pause the pro- 
gram as needed to determine intermediate values and can resume the program as desired 
at any particular point. During the program, manual adjustments may be made and run if 
desired to check certain criterion or to fine tune (step 430) for changes in the field ob- 
served values. 

It should be understood that the present invention has many advantages including 
the feature the user is able not only to calibrate using pipe roughness, but also junction 
demand and link (pipes, valves and pumps) operational status (open or closed). The pipe 
roughness can be either computed as a new roughness value for group of pipes, or modi- 
fied by multiplying the initial roughness with the correction factor. The junction demand 
can also be selected as a calibration parameter. Adjusting the junction demand enables 
the user to calibrate the model to better reflect the real demand condition when the ob- 
served data is coUected, thus it improves the accuracy of the overall model calibration. 
Further, maximum control of flie genetic algorithm optimization run time status is pro- 
vided in the present invention. As noted, a user can terminate, pause and resume the 
calibration process. 

The sensitivity of the calibrated parameter values can be investigated using a 
manual adjustment. Sensitivity analysis is the process by which model input parameters 
are varied over a reasonable range (a range of uncertainty in the value of the model pa- 
rameters) and the relative change in the model response thereto is observed. For exam- 
ple, the observed change in hydraulic flow rate or pressures can be noted. The purpose of 
a sensitivity analysis is to demonstrate the sensitivity of the model simulations to uncer- 
tainty and values of model input data. The sensitivity of one model parameter versus 
other parameters can also be demonstrated. Sensitivity analyses are also beneficial in 
determining the direction of future data collection activities. Data for which the model is 
relatively sensitive might require fiirther characterization, or more fi-equent sampling or a 
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greater number of samples over the entire network. This is opposed to the data for which 
the model is relatively insensitive, in which case the data would not require fiirther field 
characterization. 

A model can be verified to match historical field conditions. This is partly be- 
cause the choice of parameter values and boundary conditions might not be unique, in 
that more than one combination of parameter values or boundary conditions can give rise 
to the same or similar model results. History matching uses the calibrated model to re- 
produce historical field conditions. The process model of verification may result in fiir- 
ther calibration refinement of the model. After the model has successfiilly reproduced 
measured changes in field conditions, it is ready for predictive simulations. 

It should be understood that the system and method of the present invention pro- 
vides a genetic algorititei-based software process that produces calibration solutions for a 
network that is more accurate because it takes a greater number of input parameters into 
account and the method includes the ability to introduce true steady-state information into 
the calibration, and it is more efficient because it does not require repetitive trial and error 
steps. 

It should be fiirther understood, that the present invention may also be readily 
adapted to incorporate otiier parameters, or boundary conditions, not specifically men- 
tioned herein, but that may be determined to be usefiil in describing a water distribution 
network. In such a case, the present invention may be adapted to take those parameters 
into account in determining optimized calibration for the network model. 

The foregoing description has been directed to specific embodiments of the in- 
vention. It will be apparent, however, that other variations and modifications may be 
made to tiie described embodiments, with the attainment of some or all of the advantages 
of such. Therefore, it is the object of the appended claims to cover all such variations and 
modifications as come within tiie true spirit and scope of the mvention. 

What is claimed is: 
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